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Abstract 

The microscopic mechanism leading to stabilization of cubic and tetragonal 
forms of zirconia (Zr02) is analyzed by means of a self-consistent tight-binding 
model. Using this model, energies and structures of zirconia containing 
different vacancy concentrations are calculated, equivalent in concentration 
to the charge compensating vacancies associated with dissolved yttria (Y2O3) 
in the tetragonal and cubic phase helds (3.2 and 14.4% mol respectively). The 
model is shown to predict the large relaxations around an oxygen vacancy, 
and the clustering of vacancies along the (111) directions, in good agreement 
with experiments and first principles calculations. The vacancies alone are 
shown to explain the stabilization of cubic zirconia, and the mechanism is 
analyzed. 

I. INTRODUCTION 

Stabilized zirconias are materials of outstanding technological importance@. Pure zirconia 
(Zr0 2 ) is monoclinic (m) at zero temperature (space group P2i/c, Ref. 0,0); upon increase 
of the temperature (at zero pressure) the material transforms to tetragonal (t) and then 
to a cubic (c) fluorite structure (space groups P4 2 /nmc and Fm3m, respectively). These 
phase transitions induce large volume changes and make the pure material unsuitable for 
applications. The addition of lower- valence oxides like CaO, MgO, or Y2O3 disfavour the 
m phase, stabilizing more symmetric structures with cubic and tetragonal symmetryi On 
increasing dopant concentration, the material transforms to a tetragonal (if:*) form, called 
partially stabilized, and then to a cubic (c*) one, called fully stabilized. The amount of 
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doping needed for stabilization is quite substantial (no less than 8 % mol Y2O3 to achieve full 
stabilizationi^S) , and the electrostatic neutrality of this rather ionic material is maintained 
by oxygen vacancies. In this case one vacancy is required for every two yttrium ions. 

The simultaneous presence of dopant cations and oxygen vacancies in large concentration 
means that the local atomic environments in the stabilized material are very different from 
the corresponding stoichiometric (t and c) phases. Despite the analogy in the sequences 
m—t—c and m-t*-c*, there is no clear picture of the microscopic mechanisms of stabilization 
to parallel our understanding of the pure material: the most relevant issue concerns the 
respective roles of impurity cations and of oxygen vacancies^. 

We address this issue by means of a self-consistent tight-binding model, which was 
parameterized on the electronic and structural properties of pure zirconia@0. The method 
was used to calculate the temperature evolution of the free energy surfaces governing 
the c — t phase transition, which was then described within the Landau theory of phase 
transformation!! These resnlts are the starting point for the present approach to the 
stabilization mechanism, which is based on the qualitative analogy between isothermal and 
isoconcentration lines in the high-temperature region of the Zr02-Y2C>3 phase-diagram. At 
fixed concentration, higher temperatures destabilize the tetragonal phases t and £*, favouring 
the cubic ones c and c*. Similarly, at fixed temperature, higher amount of impurities stabilize 
the cubic phases. 

The position of the oxygen vacancy which is associated with each pair of Y atoms has 
been a subject of some discussion. Even though there are experimental data supporting 



the case in which Y is nearest neighbor (NN) to the vacancylBila, the most recent 
experiments0~0'0 suggest that Y is next nearest neighbor (NNN) to the vacancy. First- 
principles calculations@ii agree with the latter result, showing that the Y NNN position 
is energetically favoured by 0.34 eV with respect to the NN position.^ The analysis of 
Li and coworkers&B suggests that the Y atoms have a composition-independent 8 fold 
coordination shell (like in the perfect fluorite structure c) in both the c* and t* lattice types, 
further supporting the hypothesis that the dopants are on average NNN to the vacancies. 
The presence of oxygen vacancies, together with the Y in the NNN position, reduces the 
average coordination number of the zirconium atoms from 8, as in the c structure, to values 
closer to 7, similar to the m phase. Besides the fact that the Y atoms are 8-fold coordinated 
in both the stabilized structures, on the basis of the K edge in x-ray-absorption spectra, 
Li and coworkersEHIzI conclude that the perturbation in the neighborhood of the dopant 
cations is small and short-ranged. 

From these considerations, it is tempting to propose a physical picture in which the 
dopant cations, located in fluorite-like cation lattice sites NNN to the vacancies, do not take 
an active part in the stabilization mechanism, which is dominated by the crystal distortions 
around the oxygen vacancies. We now investigate the consequences of this hypothesis, by 
neglecting altogether the presence of the dopant atoms, therefore isolating the role played 
by oxygen vacancies in the stabilization of the c* and t* forms of zirconia. 

In Section O, the model is shown to reproduce the local atomic and electronic structure 
around an isolated oxygen vacancy predicted by first-principles calculations. The atomistic 
structures resulting from modelling partially and fully stabilized zirconia are then presented 
in Secion |TTT| , and the results are discussed in Section |IV|. Final remarks and a summary are 
made in Section |V[ 
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II. THE ISOLATED VACANCY: STRUCTURAL AND ELECTRONIC 

PROPERTIES 



As a preliminary step towards the modelling of stabilized zirconia, the structural and 
electronic properties of a model crystal structure containing an isolated oxygen vacancy 
are first addressed. The same system was studied from first principles by Stapper and 
coworkers^ a cubic 96-sites super-cell of 95 atoms in the fluorite structure containing one 
vacant oxygen lattice site, that we define as Vi. The vacancy is in the +2 charge state (one 
O 2- ion missing). In the real material, the charge of the V Q defect is compensated by the 
dopant substitutional cations Y' Zt . In our analysis, we mimic the presence of the charge- 
compensating Y atoms by distributing their charge over the 32 Zr cations, so that the cell 
is electrostatically neutral. The electronic structure in the neighborhood of the vacancy is 
described by the dangling orbitals of its neighbouring atoms. 

We define as unrelaxed the configuration in which all the atoms are in the centro- 
symmetric fluorite positions, and as relaxed the crystal structure resulting from the static 
optimization of the atomic coordinates. 

As the crystal is allowed to relax the Zr atoms NN to the vacancy radially shift outward 
along (111) directions, and the first shell of anions around the vacancy (NN O atoms) 
contract inward along (100) directions. This relaxation pattern is driven by the electrostatics 
and is shown in Figure [I]. The atomic displacements predicted by the SC-TB are compared 
to the first-principles resultsil in Table |I[ It may be noticed how the first-principles 
structural relaxations are well reproduced by the SC-TB model, without the inclusion of 
extra adjustable parameters. 

Figure ^| shows the total density of states (DOS) calculated with the SC-TB model for 
both the relaxed and unrelaxed configurations. The DOS are plotted with respect to the 
energy of the highest occupied level. Because of the choice of the reference energy, the valence 
2p oxygen bands have negative energies. Above the bandgap, it is possible to identify the 
crystal field splitting of the zirconium d bands. In the undistorted case (left panel), the 
arrow points to a sharp peak within the bandgap, also predicted by other calculations-SEI 
The corresponding partial DOS shows that this peak corresponds to the NN Zr d bands 
with symmetry E. The atomic relaxation (right panel) shifts the isolated peak back into 
the conduction band, as predicted by first principles calculations.il 

These results of the SC-TB model may be interpreted by considering the perturbation 
in the periodic electrostatic (Madelung) potential in the stoichiometric c structure caused 
by the presence of an oxygen vacancy Vq. This defect increases the Madelung potential 
on the neighbouring sites, thereby lowering the energy of the d states of the Zr ions NN 
to the vacancy with respect to the d states of the bulk Zr ions far from the defect. This 
perturbation in the Madelung potential at the nearest neighbour Zr sites is responsible for 
the isolated peak in the bandgap of the unrelaxed defective structure, marked by the arrow 
in Figure 2. The atomic relaxation effectively restores the homogeneity of the Madelung 
potential in the neighbourhood of the vacancy, pushing the NN Zr d levels back upwards. 
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III. STABILISED ZIRCONIA: STRUCTURAL STABILITY 



The study of the stabilization mechanism is based on static and dynamic simulations of 
two 96-site super-cells which were chosen as representative systems. The first one, denoted 
as Vi, contains one vacancy. This system corresponds to 3.2% mol Y2O3 and therefore it 
should be in the field of stability of the t* phase. Similarly, the second super-cell V4 contains 
4 vacancies and corresponds to 14.4% mol Y2O3, which is in the phase-field of fully stabilized 
zirconia c*. The electrostatic neutrality of the super-cells was ensured by distributing the 
compensating ionic charge of the dopant atoms over the 32 Zr cations. For each system, two 
separate sets of calculations were carried out: static relaxations and molecular dynamics 
(MD). The dynamics of the systems were simulated for 25 ps with a time-step of 5 fs, 
constraining the temperature at 300 K with a Nose-Hoover thermostatil~il. 



A. Partially stabilized zirconia 

The static minimization of the Vi super-cell has already been described. Further 
analysis of the structure reveals that the atomic configuration obtained in Section [TI| is 
metastable, corresponding to a saddle point in the potential energy. By perturbing the 
relaxed configuration obtained in Section [TT], a secondary structural modification further 
lowers the total energy of the cell, similar to the tetragonal distortion of the oxygen sublattice 
in pure stoichiometric c zirconia (soft X2 mode of vibration). 

According to the energy surface governing the c — t relationship ,@ the internal tetragonal 
distortion of the oxygen sublattice drives the external tetragonal distortion of the unit-cell. 
As a consequence, the relaxed configuration is expected to be further minimized by a non 
unitary c/a ratio. Indeed, allowing the tetragonality of the cell to adjust, the minimum 
energy is obtained for c/a = 1.02. The corresponding equilibrium structure is shown in 
Figure || The anions and cations are represented with light and dark circles respectively 
and the arrows point to the oxygen column containing the vacancy, which is behind the first 
visible oxygen ion. The projection on the x — y plane (left) shows the same pattern obtained 
with the static minimization. The NN cations move away from the vacancy and the NN 
anions contract inward. The right panel shows the projection along the x — z plane: the 
tetragonal distortion of the oxygen sublattice is very clear. Note the perturbation in the 
neighborhood of the vacancy arising from the superposition of the tetragonal distortion and 
the vacancy-induced radial displacement field. 

The MD simulation of the Vi super-cell shows that this atomic structure is stable and 
does not undergo any further structural modification. The atomic configuration obtained 
by averaging the atomic positions over the last 20 ps coincide with the one obtained with 
the static minimization, shown in Figure |3|. 

The similarity between the t* and t structures is also evident from the Radial 
Distribution Function (RDF) shown in Figure [|. In agreement with the experimental 
measurement s@'0'&"0, the splitting of the first peak in the Zr-0 RDF, a sign of the oxygen 
sublattice being tetragonally distorted, is also present in the RDF of the t* phase. The first 
peak in the 0-0 RDF is exactly the same in both the tetragonal phases, while the disorder 
introduced by the vacancy slightly modifies the second 0-0 coordination shell. The Zr-Zr 
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RDF indicates that the cation sublattice is nearly fee. Arrows in Figure |4] indicate the 
available experimental values for 3% and 15% mol Y2O3IHHEI. 

B. Fully stabilised zirconia 

The c* phase was modelled by distributing four vacancies in the super-cell. Here we 
consider two possible high-symmetry vacancy configurations. In the first one, the vacancies 
were distributed in the cell by maximizing the defect-defect separation, by setting the four 
vacancies on a fee lattice, with a defect-defect distance of \/2clq (a =9.61 a.u.=5.09 A is the 
lattice constant of the cubic structure). In the second configuration, the four vacancies were 
aligned along the (111) direction, with a defect-defect minimum distance of \f?>/2aQ. The 
SC-TB model shows that the second configuration is more stable than the first one by 0.1 
Ry per super-cell. It therefore predicts the clustering of vacancies alone the (111) direction, 
as observed experiment allyEJ'E3'E3 and from first-principles calculations." The results for the 
most stable of the two vacancy configurations, which we refer to as the V4 super-cell, are 
described in what follows. 

The static relaxation was started with all the atoms in the fluoritic centro-symmetric 
positions. In this case, only one structural modification was observed. The system 
evolved towards the atomic configuration shown in Figure [5]. The local distortions around 
the vacancies lock the c* crystal in a distorted configuration, suppressing the tetragonal 
distortion of the oxygen sublattice. The same structure was then obtained by MD simulation 
of the V4 super-cell. The local atomic environment is very different from that of the fluorite 
structure: 14 cations out of 32 are 7 fold coordinated and the majority of the oxygen ions 
are not in centro-symmetric positions. The minimum energy is obtained for a perfect cubic 
super-cell (c/a = 1). 

Despite these distortions, the RDF shown in Figure f| suggests that the structure of the 
V4 is, on average, cubic. The first peak of the Zr-0 component confirms that the oxygen 
sublattice is not tetragonally distorted. The calculated average Zr-0 distance compares very 
well with the experimental one in 15% mol yttria stabilized zirconia§Hi3, shown by the arrow. 
In the tetragonal structures, the 0-0 second coordination shell shows a characteristic three- 
peaked RDF. The corresponding RDF of the c* phase, shows only a single peak centered at 
r = a \f{2) = 6.8 a.u., which is the second shell 0-0 distance in the perfect cubic lattice. 
Because of the outward displacement of the large number of Zr atoms NN to the vacancies, 
the first Zr-Zr peak is shifted to smaller r, which is also experimentally observed (see arrow 
in Figure ^c). Apart from this detail, the Zr-Zr RDF shows complete similarity with the 
other cases and therefore the Zr sublattice is fee, even in the V4 supercell. 

IV. DISCUSSION 

Our results corroborate some of the early experimental observations on stabilized zirconia 
and on other compounds with the general formula MyO^SS These reported a strong 
distortion of the oxygen sublattice, with the NN O atoms relaxing towards the vacancy by 
« 0.3 A and the NN Zr relaxing away from the vacancy by less than 0.2 A.U Moreover, 
in the special case of the Zr 3 M 4 0i2 compounds, the difference between the ordered and 
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disordered fluorite structures was determined by the arrangement of the oxygen vacancies,^] 
which where reported to lie along the [111] directions.il 

When the results described in the previous Sections are considered in the context of the 
c «-> t phase transition in pure zirconia,0 they also provide insight on how the structural 
distortions around the oxygen vacancies affect the relative stability of the tetragonal and 
cubic phases. The relationship between the c and t forms is governed by a temperature- 
dependent double-well in the potential energy. At low temperature, the cubic phase c is 
structurally unstable and it lowers the energy by distorting the oxygen sublattice along the 
X2 mode of vibration. The extent of this distortion is measured by the order parameter S. 
The K double-well for the c — t structures is shown in Figure [| by a solid line labelled t. 

The static relaxations showed that both the unrelaxed cells Vi and V4 are unstable 
with respect to the short-range isotropic relaxation around the vacancies. Moreover, our 
atomistic simulations suggest that the Vi cell is also unstable with respect to the tetragonal 
distortion. Let us now simply superimpose these two displacement fields to see how the 
oxygen vacancies, considered as centers of dilations, modifies the energy double-well. In 
the following, the statically relaxed configurations of Vi and V4 which show the radial 
displacements only, are taken as reference systems, to which we apply the tetragonal 
distortion 5. The results are shown in Figure 

Similarly to the c structure, the defective Vi cell is unstable with respect to the tetragonal 
distortion: the atomic relaxations around the vacancy partially reduce the height and width 
of the double well (Figure ^). A static minimization from the minimum of the double- well 
leads to the same structure obtained by averaging over the MD configurations (Figure 
and showing that the finite-temperature equilibrium configuration is not just a superposition 
of the two displacement fields. 

On the other hand, the reference V4 cell is stable with respect to the tetragonal distortion 
of the oxygen sublattice and there is no double- well in the total energy (Figure §). The 
energy profile is highly anharmonic, just as the high-temperature one in the free-energy for 
the pure phaseH. This means that in our model for c*, it is possible to shift entire column of 
oxygen atoms by 5-10 % of the Zr-0 interatomic distance at no energy cost, with important 
implications for the mobility of defects and therefore for transport phenomena in stabilized 
zirconia. We know from our previous simulations, and it is not surprising, that in a flat well 
of this form the average structure at finite temperature is cubic. 



V. CONCLUSIONS 

In conclusion, partially and fully stabilized zirconia were modelled by introducing 
respectively one and four vacancies in 96-site cubic super-cell (Vi and V4, respectively). 
Energetics and atomic structures of the defective super-cells were calculated with a SC-TB 
model which was parameterized on the properties of pure stoichiometric zirconia. 

The analysis of the Vi and V 4 equilibrium crystal structures suggests a possible 
explanation of the stabilization mechanism. When the concentration of vacancies is low, 
as in the Vi super-cell, a relatively large volume of crystal is left in the fluorite structure, 
and it undergoes the tetragonal distortion just as the stoichiometric material in the c-phase 
would. Because this distortion involves the coordinated motion of all the oxygens, also the 
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atoms neighboring the defects are dragged along the tetragonal distortion of the oxygen 
sublattice (see Figure 

When the concentration of defects is high, as in the V4 super-cell, there is virtually no 
undistorted cubic region in the statically relaxed super-cell; every oxygen atom is either itself 
a neighbor of a vacant site, or at least four of its six neighboring oxygen atoms are. Thus 
there is virtually no region in which the local atomic environment is close to fluorite, which 
could undergo the tetragonal distortion (see Figure ^), and the radial distortions around the 
vacancies dominate the finite-temperature equilibrium configuration. The resulting atomic 
structure is cubic only by averaging over relatively large number of atoms (~ 100), but the 
short-range atomic structure does not have cubic symmetry. 

These results show that the stabilization of the t* and c* structures may be achieved in 
theory by doping zirconia crystals with oxygen vacancies only, and support the idea that 
the electronic and structural properties of stabilized zirconia are controlled by the structural 
disorder around the oxygen vacancies, rather than by the substitutional dopant cations. 
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FIGURE CAPTIONS 

FIG. 1. Crystal relaxation of the atoms neighbouring an isolated vacancy in a fluorite 
lattice. Key: Zr black, O dark grey vacancy light grey. 

FIG. 2. Total density of states for the relaxed and unrelaxed configurations of the V\ 
supercell, calculated with the SC-TB. The of energy is the highest occupied electron state. 
The arrow points to the energy levels of the Zr atoms NN to the vacancy split from the d 
bands of the "bulk" Zr atoms. 

FIG. 3. Equilibrium crystal structure of the Vi supercell. The arrow point to the oxygen 
columns containing the vacancy. 

FIG. 4. Equilibrium crystal structure of the V 4 supercell. The arrows point to the oxygen 
columns containing the vacancies. 

FIG. 5. Radial Distribution Function g(r) obtained from the MD simulations at 300 K for 
three supercells: t-Zr0 2 , Vi and V 4 . Arrows represent the experimental values for t*- and 
c*-Zr0 2 (Ref. 13). 

FIG. 6. Total energy vs. tetragonal distortion S of the oxygen sublattice for three 96-site 
cells: t pure stoichiometric tetragonal structure; t* partially stabilized zirconia (Vi super- 
cell); c* fully stabilized zirconia (V4 super-cell). 
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TABLES 

TABLE I. Atomic displacement (in A) predicted by the SC-TB model compared to ab initio 
results. 1 





NN 


Zr NN 


NNN 


Zr NNN 


DFT 
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Ref. 
this 1 


28 

NO 


rk 


0.24 
0.27 


0.14 
0.18 


0.04 
0.03 


0.04 
0.03 



x The displacements are along the directions shown by the arrows in Figure 1. NN and 
NNN mean Nearest Neighbour and Next Nearest Neighbour to the defect. 
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FIGURES 



FIG. 1. Crystal relaxation of the atoms neighbouring an isolated vacancy in a fluorite lattice. 
Key: Zr black, O dark grey, vacancy light grey. 




FIG. 2. Total density of states for the relaxed and unrelaxed configurations of the V\ supercell, 
calculated with the SC-TB. The of energy is the highest occupied electron state. The arrow 
points to the energy levels of the Zr atoms NN to the vacancy, split from the d bands of the "bulk" 
Zr atoms. 
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FIG. 3. Equilibrium crystal structure of the Vi supercell. The arrow point to the oxygen 
columns containing the vacancy. 




FIG. 4. Radial Distribution Function g(r) obtained from the MD simulations at 300 K for three 
supercells: t-Zr02, Vi and V4. Arrows represent the experimental values for t*- and c*-Zr02 (Ref. 
13). 
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FIG. 5. Equilibrium crystal structure of the V4 supercell. The arrows point to the oxygen 
columns containing the vacancies. 
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FIG. 6. Total energy vs. tetragonal distortion 8 of the oxygen sublattice for three 96-site cells: 
t pure stoichiometric tetragonal structure; t* partially stabilized zirconia (Vi super-cell); c* fully 
stabilized zirconia (V4 super-cell). 
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